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^\ , Recent developments in the analysis of finite temperature dissipationless transport in inte- 

grable quantum many body problems are presented. In particular, we will discuss: (i) the 
f— I I role played by the conservation laws in systems as the spin f/2 Heisenberg chain and the 

r—{ . one- dimensional Hubbard model, (ii) exact results obtained using the Bethe ansatz method 

on the long time decay of current correlations. 
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1 Introduction 

It has recently being proposed that integrable quantum many body systems show dissipationless 
Tj" \ transport at finite temperatures El. This idea was motivated by analytical and numerical studies 

''^ ■ on a toy model □ and the fermionic version of the Heisenberg modelcl. It is analogous to the well 

f— ^ . known effect of transport by solitons in classical nonlinear integrable systems. In this domain of 

0\ \ classical physics, a part of the activity is already in the technological applications, for instance 

^ ■ the propagation of solitons in optical fibersQ and in particular the robustness of dissipationless 

"j^ . transport to perturbations. 

In quantum many body systems, there is not yet clear experimental evidence for unusually 
high conductivity (close to ideal) at finite temperatures for a system described by a (nearly) 
^ ■ integrable Hamiltonian. Perhaps the best candidates so far are quasi one dimensional spin 1/2 

O . systems described by the Heisenberg model. In fact NMR studies^ on the Sr2CuO^ compound 

showed an unusually high value of the diffusion constant, characteristic of ballistic rather than 
diffusive behavior. 

^ , Besides the magnetic compounds, systems characterized by strong short range electronic 

H \ correlations (described by the integrable one dimensional Hubbard model) could be candidates 

for observing the predicted ballistic behavior, even more conspicuous at high temperatures, 
instead of the expected diffusive one due to electron-electron scattering. Nanotubes or artificially 
made nanostructures come to mind as possible experimental realizations. 

From the experimental point of view, the most relevant question is the robustness of this dis- 
sipationless transport to perturbations, for instance impurities, 3d coupling, phonons. The same 
problem analyzed for classical systems, indicates relative insensitivity of the soliton propagation 
to disordercl. From the theoretical perspective, an analytical solution of the dynamic properties 
of integrable systems should be expected, exactly due to the integrability of the models. On the 
other hand, the analysis of robustness might largely depend on numerical simulation studies, 
the same situation as in classical systems. We should also point out, that this problem is closely 
related to studies on dynamical systems, analyzing the conditions of appearance and stability 
of classical/quantum chaos Q'LI. 

The framework for discussing dynamic properties at finite temperatures/frequencies is the 
linear response theory (or Kubo formalism). In this formulation, using the fluctuation - dis- 
sipation theorem, the transport properties are related to the study of dynamic correlations at 
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equilibrium. It might be objected that results drawn from this method for integrable systems are 
questionable. However, we can point out that transport pronerties derived by linear response 
theory for noninteracting systems (e.g. an harmonic crystal) eI or classical integrable systems 113 
are in accord with experimental observations u and results from more sophisticated theoretical 
methods. 

The basic way for characterizing the transport properties is by the dynamic current cor- 
relations < j(t)j >. As a criterion for dissipationless transport we will study the long time 
asymptotic value Cjj, < j{t)j{0) >t~KX]= Cjj. If Cjj is finite, the response is reactive and the 
system shows ideal conducting behavior. Cjj is related to the Drude weight D(T) in the real 
part of the conductivity as expressed in linear response theory: 

a{u) = 2TTD6iu;) + aregico), D = ^Cjj (1) 

It should be noted that a vanishing D{T) does not necessarily imply diffusive behavior. If 
the current-current correlations decay too slowly then transport coefficients cannot be defined. 
In the analogous situation of one dimensional classical systems, different behaviors have been 
observed t3; typically if a system is integrable then Cjj > and the system is an ideal conductor. 
If it is nonintegrable, then models have been found where correlations decay fast enough and 
a transport coefficient can be defined, but there are also models where the decay is too slow 
indicating anomalous low frequency conductivity. In quantum many body systems the situation 
is less clear, but we might similarly expect a variety of behaviors. 

In the following we will discuss two recent developments for analyzing the transport proper- 
ties in integrable systems: the first, it is the way in whidi the conservation laws characterizing 
these systems affect their long time current correlations Eil; the second, it is a procedurell3 based 
on the Bethe ansatz method for calculating analytically D{T). 

2 The role of conservation laws 

Integrable many body systems are characterized by a macroscopic number of conserved quan- 
tities 113. A set of conservation laws is represented by local involutive operators Qn, commuting 
with each other [Qn,Qm] = and with the Hamiltonian, [Qn,H] = 0. We can relate the time 
decay of correlations to local conserved quantities in Hamiltonian systems by using an inequality 
proposed by MazurO 

hm ^ r < A{t)A >dt>y ^ ^^" ^' (2) 

Here <> denotes thermodynamic average, the sum is over a subset of conserved quantities Qn, 
orthogonal to each other < QnQm >=< Qn > ^n,rm A^ = A and < A >= 0. Applying this 
inequality to current-current correlations we obtain: 

We will now discuss two applications of this inequality, the decay of spin currents in the 
Heisenberg model and the electrical conductivity in the Hubbard model HH. 

2.1 s=l/2 Heisenberg model 

It is described by the Hamiltonian: 

L 

H = 22{JxSfSf_^_l + JySfSf^^ + J^SfSf^i) (4) 

1=1 



(by a Jordan- Wigner transformation the Heisenberg model is equivalent to a model of spinless 
fermions interacting with a nearest-neighbor interaction (t — V model)). 

This model is characterized by a macroscopic number of conservation laws. It is worth to 
note that the first nontrivial conservation law Qs is related to a physical quantity, it coincides 
with the energy current j^: 

Qs =f = T. J.JyiSf-iSfSl, - Sl.SfSf^,) + (x, y, z). (5) 

1=1 

This is particularly interesting and of actual experimental relevanceEJ as it implies ideal thermal 
conductivity (energy currents do not decay at all, < j^{t)j^ >= const.). It should also be 
taken into account in the analysis of the nuasielastic peak in Raman experiments where the 
local dynamic energy correlations enter EJIIZI. 

Regarding the spin current correlations, for Jj = J^ = 1, J^ = A we find that: 

D = Ac,.,. > Al£Qi2l (6) 

where j^ =spin current= J2i=i{Sf Sfj^i — SfSfj^^), implying ideal spin conductivity and the 
absence of spin diffusion. We can analytically calculate the right hand size of this inequality in 
the /3 — > limit, 

^ /? 8A2m2(l/4 - m^) 

D > -7--S TT, m =< Sf > (7) 

From this result we notice that we obtain a positive bound only for a finite magnetization m, 
when the system is in a magnetic field. For 171 = this inequality does not provide a useful 
bound although, as we will show below using the Bethe ansatz method, D(T) is still positive. 



2.2 Hubbard model 

Finally, for the Hubbard model described by the Hamiltonian, 



L ^11 

H = i-t) Y. (cLq+i^ + h.c.) + t/^(n,t - -)(n,; - -) (8) 

a,i=l i=l 

considering Q^ given by: 

Q3 = "Yi-tfiicl+ia^^i-ia + h.c.) - C/(ji„i,i,^ + ji,i+i,CT)(ni _^ - -) (9) 



we obtain for the decay of current correlations in the /3 — > limit, 

r, > P [£/E<.2p.(l-pJ(2p_,-l)]2 

- 2 E.2p.(l - p^)[l + U^2pl^ - 2p_. + 1)] ^ ^ 

Again, we obtain a finite value for D(T) and so ideal conducting behavior for a system out of 
half-filling. Bethe ansatz analysis by Fujimoto and Kawakaniill3 confirm this result (the system 
at half- filling has also been analyzed). 



3 Bethe ansatz analysis for the Heisenberg model 



Recentlva new procedure was proposed for analytically calculating D{T) using the Bethe ansatz 
methodt3.The application to the Heisenberg spin 1/2 model, that we will now sketchOJ, is based 
on a series of assumptions and developments. 

First, the structure of the solutions of the Bethe ansatz equations for < A < 1 conjectured 
by Takahashi and Suzuki EJ is used. Parametrizing the anisotropy coupling by A = cos 6, 9 = 
^, u integer, and the pseudomomenta k^ by the rapidities a;cM--Pot(-2^) = cot(2) tanh(-|^), the 
solutions are organized into a finite number of allowed strings c3, n = l,2,...,i/ — 1, 



n,k 



Xn 



xl + {n + l-2k)i + 0{e 
Xa + iv + 0{e 



'SN\ 



k = 1,2, ...n 



-6N\ 



6>0. 



(11) 
(12) 



Next, finite size corrections on the positions of the strings are expressed ciJ by introducing 
the functions gij and g2j'- 



'^N 
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(13) 



where xjy(x^) are the rapidities for a system of size N{oo). 



Finally introducing excitation and hole densities pj and Pj as in the standard thermodynamic 



Bethe Ansatz and expanding the Bethe ansatz equations to 0{jj') and 0{jt^) we obtain: 



^=of^Y. dx{pj + p]) < nj > (1- < Uj >)( 



9^j dgij 
dx dd> 



|(/>^0 



(14) 



where < rij >= 1/(1 + e^'^^) and ej = {1 / P) ln{pj / pj) . This expression looks formally analogous 
to that for free fermions: 



D 



JL 

2N 



de„ 



5] < n^ > (1- < n,, >)(-£)f^_^o 



(15) 



where < n^ > is the Fermi-Dirac distribution and e^ are single particle energies. 

The functions pj, p'j, gij, g2j are determined numerically by iteration. Using this procedure 
we obtain D{T) for all temperatures and couplings < A < 1 as is shown in Fig. 1. 
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Figure 1: 0(A) evaluated at the points v ■ 



16 and various temperatures. The continuous line is the high 



temperature proportionality constant Cjj = D/f3. The o indicate exact diagonalization results from ref. (17) 



In particular, we recover the known T = value E3 Dq with a characteristic power law 
behavior: 

D{T) =Do- const.T'', a = 2/{u - 1), (16) 

and we obtain good agreement with numerical simulations O for f3 = (diamonds), a result that 
lends support both to the Bethe ansatz procedure and the string assumption. 

Particularly interesting is the vanishing of D{T) for A = 1 at any finite temperature. This 
result suggests that for A > 1 the Heisenberg model does not show ideal conducting behavior 
and the isotropic model is a borderline case. If the conjecture, born out of a numerical studya, 
of absence of weight in the conductivity at low frequencies proves valid in the thermodynamic 
limit, then we would have a realization of an exotic ideal insulating phase. Finally, from the 
results presented, it is not unreasonable to expect anomalous low frequency spin dynamics for 
the isotropic Heisenberg model, hopefully observable in experiments. 
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